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Abstract 

In the framework of scalar-vector dominance we calculate the hadronic matrix el- 



' ements of scalar and tensor effective quark currents induced by virtual leptoquark in- 

■ teractions. Combined bounds on the product of couplings and leptoquark masses are 

Oh| obtained from experimental data. 

^ ; 1 Introduction 

m _ 

^ ■ Recently, new data on the analysis of -k^I^vi decays were published by two collabo- 

\Q '. rations: KEK-E246 |l[ and ISTRA 0, in addition to the presentation given by the Particle 

O \ Data Group 0. So, at present we have got quite precise measurements of characteristics in 

Q I the Ki'i decays, that needs a theoretical interpretation in the framework of Standard Model 

CN ■ (SM) as well as beyond it. Such the study is of interest because of the experimental search 

^ for effects, which can point to the contributions with the violation of combined CP-parity in 

. the kaon decays, for example, the transverse T-odd polarization of lepton in Ki^^ modes 

I \ that can essentially enrich the information on the CP-breaking dynamics in addition to the 

g_) \ program with the B- mesons 0. 

^ • The matrix element of decay is parameterized in terms of scalar, vector and tensor form 



factors, fs, f± and /t, in the following general form 



L niK 
where the lepton currents are given by the expressions 



'11 



so that chiral spinors are 



^R.L = ^(1±75)^, 



and Gp is the Fermi constant, rriK is the mass of kaon. The four-momenta are defined as 

P = Pk+P7t, q=PK-PiT, 
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while Vsu is the matrix element of Cabibbo-Kobayashi-Maskawa matrix for the mixing of 
weak charged quark-currents. We define the generators a by the commutator 

The dependence of form factors on is usually expressed in terms of linear slopes normalized 
to get the dimensionless quantities 



^ _ dln/,(g2 



dg^/m^ 

where m.,^ is the pion mass. The combination of form factors 



(2) 



/o = /+ + ^/-, (3) 
p-q 

is introduced, so that the experimental data are given in terms of the following set 0: 

fs T 

The data of [0, H on A+,o can be averaged, so that with the statistical errors we get 

A+ = 0.0287 ±0.0018, (4) 
Ao = 0.0203 ±0.0033, (5) 

while the systematic uncertainties are given in the original papers. The values in @) and (^) 
result in the ratio 

Zi-M = _o.096 ± 0.043. (6) 
/+(0) 

The given parameter A+ is in a good agreement with the PDG values for both the electron 
and muon modes 0, while Aq and /_(0)//+(0) above are within the limits of 1.5 a-deviations 
from the PDG averages. The preliminary analysis by KTeV gives 

A+ = 0.0275 ± 0.0008, 

which is close to the estimate in (^). 

In the framework of SM we get the form factors 

(vr°(p.)| S7^« \K+{Pk)) = ^ {Up, + f-q,), (7) 

while 

fs = fT = 0. 

Therefore, the study of scalar and tensor form factors is a good test for the search of 'new' 
physics beyond the SM. 
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Supposing we derive 



/o 



{7i\p^)\ -s^,u \K+{pk)) = {mu - m,) {n\p^)\ 7su \K+ipK)) = (p ■ q) ^, (8) 

implying that the form factor /o determines the matrix element of scalar quark-current. 
The contraction of vector lepton-current 

-l^^q^f- = mif.ls, 

induces the scalar term in the matrix element. So, the electron mode is more sensitive to the 
extraction of scalar form factor, since the SM background contribution is suppressed by the 
lepton mass, 

fSM _ f 

At present, the measurements of fs and /t result in values slightly deviating from zero, that 
is consistent with the expectations of SM. So, in the electron mode 

fs 



/+(0) 
/t 



0.0040 ± 0.0160(stat.) ± 0.0067(syst.), (9) 
-0.019 ± 0.080(stat.) ± 0.038(syst.), (10) 



/+(0) 

where we have taken into account the redefinition of sign in comparison with the appropriate 
formula in |I| as accepted in this paper in (p, while the combined analysis of muon and 
electron modes in results in the similar values 

fs 



/+(0) 
/t 



0.004 ± 0.005(stat.) ± 0.005(syst.), (11) 
-0.021 ± 0.028(stat.) ± 0.014(syst.). (12) 



/+(0) 

The collaboration KTeV presented the following constraints in the electron mode 

fs 



/+(0) 
/t 



< 0.04, (13) 

In the present paper we study nonzero contribution of leptoquark interactions to the tensor 
form factor, which correlates with the scalar one due to the Fierz transformation. In section 2 
the effective lagrangians with the virtual leptoquarks are described as concerns for the decays 
of —>■ TT^l'^vi, and the requared matrix elements of quark currents are presented. General 
expressions for the hadronic matrix elements with the tensor structure are derived in section 3, 
where we develop the model based on the dominance of vector and scalar mesons and adjust it 
in the description of f±^o form factors. In the framework of potential approach the preferable 
region of model parameter is limited in agreement with the experimental data. The constraints 
on the masses of scalar leptoquark and their couplings to the fermions are obtained in section 
4. The results are summarized in the Conclusion. 
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2 The contribution of leptoquark interactions 



A consistent classification of leptoquarks under tlie gauge symmetries of SM were done by 
Buchmiiller, Riickl and Wyler in ^j. We accept the nomenclature prescribed in as shown 
in Table |l] extracted from |10|. So, the leptoquarks are marked by their spin, representation 



of weak 5't/(2)-group (singlets, doublets and triplets), appropriate electric charges in the 
multiplets and the fermion number F. For the sake of briefness, the flavor of lepton is marked 
by the electron in Table |1|, while the couplings Yl.^ should be labelled by the flavor indices, 
too. 

The diagrams describing the contribution of leptoquark interactions into the form factors 
under study are shown in Fig. |[ 




Figure 1: Two kinds of leptoquark exchanges contributing to the tensor form factor in the 
decay ti^I^ui. 

The tensor terms appear under the Fierz transformations, so that the vector leptoquarks 
do not contribute into the tensor form factor. Further, the tensor term shifts the helicity of 
leptons. Therefore, we isolate the leptoquarks involving the interaction with both the left- 
handed neutrinos and right-handed charged leptons. The appropriate vertices are shaded in 
Table |l]. Thus, we consider the following scalar leptoquarks: the singlet 5*0 and the doublet 
5*1/2 with the charge —2/3. 

The Yukawa-like interactions involving the strange quark have the form 

C[S,/2] = Sl/^{Y^u^u^ + Y^sJ^)+h.c., (15) 

C[So] = S*[Y[°](Mc,aeL + Sc,R,/^L)+Y[°UcA]+h.c., (16) 

where we have omitted the flavor indices. These lagrangians induce the effective low-energy 
interactions according to the formulae 



1 Y Y* 
o IVIq 

*l/2 

--^S^(^LnR)(z?Ja)+h.C., (17) 



2 Mi 
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Table 1: The first generation scalar (S) leptoquarks/squarks and vector (V) leptoquarks in 
the BRW model p| according to the nomenclature in P| with their electric charge in units of 
e and fermion number F = L + 3B. For each possible non-zero coupling Y the decay modes 
and the corresponding branching ratio /3e for the decay into an electron and a quark are also 
listed. The restrictions on the values of (3e arise from the assumption of chiral couplings. 
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where we have used the Fierz transformations for the chiral fermions, taking into account the 
identity 

75 = ^'r, 

that causes the summation of scalar and pseudoscalar parts (the factor of 2). The anti- 
commutation of fermions has been explored, too (the overall negative sign). Further we 
introduce the notation 



since the above combination of leptoquark mass and couplings enters the problem under study. 

As for the contribution of 5*0, one can easily find that the effective lagrangian has the same 
form of (|17|), because the charge conjugation of spinor is defined by 



9c = ce*, 

where C = 172 in the Dirac representation of 7-matrices, so that 

where T denotes the transposition. The terms induced by the leptoquarks 5*1/2 and 5*0 can 
interfere, of course. However, we include this effect into the definition of scale A^q. 

Thus, we can estimate the contribution of leptoquark interactions, once we calculate the 
appropriate matrix elements of quark currents, that is the deal of next section. 



3 Hadronic matrix elements 

The experimental data on the slopes of form factors shown in the Introduction are in a good 
agreement with the estimates in the framework of chiral perturbation theory (xPT) ||1 1|| . 
However, to the moment we have not any predictions of xPT on hands as concerns for the 
hadronic matrix elements of tensor quark-current. In the present paper we explore the model of 
meson dominance, i.e. the dominance of vector and scalar states appropriate for the quantum 
numbers of transitions between the quarks. The corresponding diagram is shown in Fig. |^. 
Considering the vector quark-current, we can evaluate the form factors 

r I 2\ Ik* rriK* 1 . x 

/-(?^) = -/+(?^)^%^ + ^x*x.4^^ (20) 



in terms of couplings entering the following LagrangiansQ 



f^K*K-K = gK'KTviPK + P7t),i(^k* Vk^Ij (21) 
^K*K7r = gK*KnVKSVKVl^ (22) 



^The couplings g are prescribed for the charged pions, while the neutral ones have the isospin factor l/\/2- 
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TT 



Figure 2: The diagram describing the contribution of excited vector and scalar kaon states 
into the hadronic matrix element of current j factorized from the lepton part in the decay 

and 

{K*{k)\-s^^u\0) = fK^efniK*, (23) 
{K;{k)\-s^,u\Q) = fK^^k,, (24) 

where ip denotes the appropriate field, and e is the polarization vector of K*. In (|19]) we 
have introduced the running pole mass msuiff) in the transition s ^ u. The normalization 
condition is rather evident 

while we need the value at = 0, rrisu = rnsu{0), since we use the approximation of linear 
evolution of form factors, 

U{q') ^ 9K*Kn^^^^{l + qVml), (25) 
f-il') - -Uici')^^^^^^ + 9K*K.^il + qVrnl.). (26) 

The evolution of msu{<f) to = is expected to be slow in the framework of model with the 
meson dominance. We suppose that the spin forces in the bound state should be suppressed 
beyond the pole, since they depend on a density of bound states, which drops outside the 
poles. So, the spin-averaged mass of IS-level in the su system is known experimentally, 

msu[lS] = ^{rriK + SrriK*) ~ 793 MeV. 

We expect that 

[IS] < rUsuiO) < rriK*. 
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So, we put 

rrisu ~ \{m,u[lS] + m^.) ~ 0.85 GeV, (27) 

which is inside the systematics uncertainty of the model. Since the spin-dependent forces are 
suppressed in the excited P-waves, we put the pole mass in the scalar sector to be equal to 
the experimental value of Kq. 

From (|2^) and (^) one can easily deduce the expression for the scalar- channel form factor, 

/o(g^) - /+(0) + (28) 



as well as the slopes. 



where 



A+ = ^, (29) 
Ao = 5-A+, (30) 



gK^KnfK^ml^ ^^^^ 



/+(o) ^Ksi^K -^D' 

/+(0) = gK^K.^^^^. (32) 

The most of model parameters can be extracted from the experimental data. So, the coupling 
constant fx* is well known, 

^ 215 MeV, 

while the decay constants g are related with the widths measured^. 



r[R-^K^] = A-K,^^, (33) 
TK^K.] = (34) 

whereas \Pk\ denotes the momentum of kaon in the c.m.s, so that numerically^ 

gK*Kn ~ 3.94, gK*Kn ~ 3.48 GeV. 

^In the formulae for the total widths of K* and Kq we have explored the isospin-symmetry relations: 
T[K*+ K+n^] = l/2r[K*+ -> K°t:+] and the similar equation for the scalar meson K^. 

•^In the estimates we put the effective masses in the equations relating the constants with the total widths, so 
that mx" — > rrisui^S] and mx* — 2|p^| in the limit oirriK* 3> rriK, tttLt^. In the phenomenological model under 
study, the decay constants g enter the form factors in terms of products with the leptonic couplings /. These 
products should be adjusted in order to satisfy some conditions motivated by QCD and its chiral symmetry. 
In this way we have to follow a specified approach in estimates for both g and / as given below. We stress that 
the model parameters g are quite uncertain because of reasons inherent for the phenomenological approach 
ignoring higher excitations as well as a continuum contribution. Neverteless, we argue for the preerable choice 
of numerical values. 



|3 
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The only free parameter of the model is the coupling fx*, which we are tending to restrict in 
the framework of potential calculations by the comparative analysis with the known leptonic 
constants of p and K*. For this purpose, we calculate the diagram in Fig. ^, where the 
quark-meson vertex includes the wave function of constituent quarks. 



K*, K* 




Figure 3: The diagram describing the contribution of quark loop into the hadronic matrix 
element of current j. 

For the scalar state we use the current 
with the identity 

idf, [s{x)Yu{x)] = {mu - ms)j{x). 

In this technique we find 

f!q = ^^^^:^z-^L=I^L(0)|, (35) 



/™ = \ r7^,\RU0)\, (36) 

where -Rsu(r) denotes the radial wave function in the system of su, msul^S] is the spin-averaged 
mass of K* and K, and mred is the constituent reduced mass for K^, so that 

« '"t;'^f!""-'\fL « 0.34 GeV. 



For the p meson we have the expression similar to (pq) under the substitution s ^ d. 

Further, we explore the static potential derived in |jl2| and solve the Schrodinger equation 

— + V{r) 



^{r) = [A{p,) + 2{po-fi,)]^{r), (37) 



for the system du, so that the binding energy A{pq) of IS'-level is related with the mass 

msu[lS]=A{p,)+26p, (38) 

and it is shown in Fig. Ijat /io = 0.345 GeV, /i* = 0.224 GeV versus the light quark constituent 
mass pq with 6p = p* — po- In (^8]) we do not add the constituent masses of light quarks into 



mdui^S], GeV 
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Figure 4: The mass of 15* level in the system du calculated in the potential model with the 
constituent mass /iq. 



the mass of meson, since the constituent masses are really the parts of potential energy V{r) 
in the confining quark-gluon string. 

The mass of bound state shows the minimum versus the constituent mass at /i*, which 
gives the optimal value of mass for the calculation of radial wave function. For the constituent 
mass of strange quark we use 

Hs=ms+ Hg, 

with rris = 0.24 GeV, which represents the current mass at the scale of 1 GeV |T^. 

At this stage the estimates of coupling constants in the potential model can be optimally 
got according to ( P^D and (^). However, the corrections by both the quark-gluon loops and a 
relativistic motion can be rather essential, that can be taken into account by the introduction 
of /C-f actor, 



At 



we get the estimates 



1.45 



fp 
Ik* 



205 MeV, 
217 MeV, 
130 MeV. 



(39) 
(40) 
(41) 



The /C-factor should generally depend on the spin and flavor of current under srudy. The above 
estimates show that the dependence on the flavors of quarks composing the bound state is 
rather suppressed, since we have amazingly reproduced the coupling constants of vector states 
in the limits of experimental intervals with the uniform /C-factor. As for the dependence on 
the quantum numbers of the meson, we expect that the variation of /C-factor is negligibly 
small because the summed quark spin in both K* and Kq is equal to 1, while the spin-orbital 
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contributions are usually suppressed. Thus, the estimate in (l4l|) should be quite accurate up 
to 5 MeV, as it does for p and K*. Nevertheless, we permit a conservative variation^ 



120 MeV < fK'^ < 140 MeV. (42) 

Further, we can compare the model estimates with the experimental data on Ki^ decays 
listed in the Introduction. This analysis is presented in Figs. ^ and |^. We draw the conclusion 
on the model is well adjusted in describing the data. 

Ao 

0.03 
0.02 




0.01 



0.1 



0.12 



0.14 



0.16 



GeV 



Figure 5: The model predictions for the slope Aq versus the coupling constant of meson 
(the solid line) in comparison with the experimental data (the horizontal band). The vertical 
band gives the region of preferable values of fx* expected from the potential model. 



According to ( [29| ) and ( p2D the values of A+ and /+(0) are independent of fx*- Numerically, 

we get 

A+ = 0.0271 ±0.0011, /+(0) = 1.046 ±0.040, (43) 

which are in a good agreement with both the experimental data and predictions of xPT. 

Further, we test the model under the Callan-Treiman relation, that expresses the sum of 
vector-current form factors in terms of leptonic constants of kaon and pion: 

/+«)+/-«) = ^- (44) 

In the model under study we get 



See the sum rule estimates in [13 



11 




Figure 6: The model predictions for the ratio /_(0)//+(0) versus the couphng constant of Kq 
meson (the sohd hne) in comparison with the experimental data (the horizontal band). The 
vertical band gives the region of preferable values of fx* expected from the potential model. 

where we have neglected the kaon mass with respect to the scalar meson one. Then, numeri- 
cally the relations result in 

^ = 1.305 ±0.020 or ^ ^ 1.273 ± 0.017 

/tt Jit 

under the variation in (^2|). So, at = 132 MeV we deduce 

fx = 172 ± 3 MeV or fx = 168 ± 4 MeV. 

The systematic error caused by the approximation, as we see, is about 5 MeV, and conserva- 
tively one expects 

/ic = 170 ± 4 ± 5 MeV, 

which is in a good agreement with the known data. 

Neglecting both the deviation of /+(0) from the unit and the kaon mass with respect to 
the mass of scalar Kq , we can derive from (^4D and (|45|) the Dashen-Weinstein relation 

m\-ml ( ^) ■ 

Both relations by Callan-Treiman and Dashen-Weinstein can acquire valuable numerical cor- 
rections in the model under study as well as in the xPT. So, from the formula for Aq we 
get 

fK = 160 ± 4 MeV, 
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so that the displacement of fx points to the possible size of corrections. 

Then, we calculate the expression for the hadronic matrix element of tensor quark-current 

(7r°(p^)| sa^^u \K^{pk)) = -i -^^^{p^^Qu - PuQf,) ~ -i ^^^(p^g^ - Puq^,), (46) 

where we have neglected the dependence of /+ on g^, since the antisymmetric tensor is linear 
in q. Formula (^) can be compared with the general expression 

(7r°(p^)| sa^^u \K+{pk)) = -^B {p^q^ - p^q^), (47) 

so that B depends on a single additional quantity c_(g^) 

B = C- h [ms + nriu) 



ms-rriu p-q 
where we have explored the definition 



(7r°(p^)| i|s(9^M) - {d^s)u^ |ir+(j9i^)) = {c+p^ + c_ g^) (7r°(j9^)| su \K'^{pk)), 



with an evident condition of self-consistency 

c+ = (m, -ni^ + c-q ). 

p-q 

In the model of meson dominance we get 

^ _ f+ rris- rriu /_ - ml 
fo rnK" fo p-q 

Neglecting the current mass of light quark, at = we find 



(4^ 



c-(0) ^ ^ + (A+-Ao)4' (49) 

c+(0) ^ 2^^. (50) 

The physical meaning of c± is rather simple: they determine the difference between the frac- 
tions of meson momenta carried by the s and u quarks in the kaon and pion under the weak 
transition. At g^ = we get 

aK = ^(c+ + c_) ^0.018, (51) 

= ^(c_-c+)^ 0.28. (52) 
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Spin effects: the polarization of lepton in SM 

Neglecting the suppresed contributions by the scalar and tensor form factors in (|l]), we cal- 
culate the matrix element squared for the lepton with the spin polarization Sa, satisfying the 
conditions of = — 1 and s ■ pi = 0. Then, omiting an irrelevant overall normalization factor 
and putting the vector form factors to be real, we get 

1^ IK+ ^ TT^Z+z/,] p ~ 2fl {pi ■ p) {p, ■p) + 2Uf^ {p, ■ p) + {pi ■ p,) {f_ mf - fl p') - 

'^f+f- m (pi ■ Pu) (p- s)- 2fl mi{p^ ■p)ip- s)+ (53) 
2/+/- m {pi ■ p) {p^ ■ s) + mi fl + f_ mf) {py ■ s). 

Following the ordinary definition for the polarization Pj, 

\ i=L,N,T 

where denotes the longitudinal four- vector 

ex = ^--{\Pi\\EiPi), 
m\Pi\ 

the unit vector lies in the decay plane, and it is orthogonal to e^, while is transversal 
to both cl and e^v. 

The scalar products in (|53|) can be easily expressed in terms of pion and lepton energies. 
So, for the longitudinal polarization we get 

Pu-s = [Ey\Pi\ - EfEy + El {pi ■ p^)], 

mi\Pi\ 

P-s = -^[EM-EfE^ + Ei{pi-p^y ' ""^'^'l 



mi\Pi\ mi 

The results for the muon and electron modes are shown in Figs. ^ and ^, respectively. 

As we can expect from (^3]), the variation of longitudinal polarization is significant for the 
muon in contrast to the electron, since the spin dependence exhibits the suppression by the 
lepton mass. The longitudinal polarization effects for the electron mode are negligibly small. 

In Fig. we see that the longitudinal polarization of the muon is essentially negative in 
the dominant part of Dalitz plot. The significant variation of form factor /_ by a factor of 5 
shows a weak sensetivity of longitudinal polarization to /_. Fig. | exhibits that the deviation 
of electron polarization from —1 is negligible. 

Further, we have calculated the normal polarization of leptons as shown in Fig. |^. We see 
that the muon has a significant variation of normal polarization, while the electron spin in 
the normal direction is quite small (about 0.2%), and it is sizable (about 10%) in the region, 
where the number of events is essentailly suppressed. 
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0.5 0.6 0.7 0.8 0.9 



y = 2Ef,/mK 

Figure 7: The distribution of longitudinal polarization on the Dalitz plot for the decay of 
—>■ Ti^fi'^h'^ with X = 2ET,/mK, y = 2E^/mK (the left picture). The contours are shown 
with the arithmetic progression for the polarization: —0.9375 +j ■ 0.125. The section of Dalitz 
plot a,ty = 0.8 with the model value of /_ described in the text (the solid curve) in comparison 
with the fixed /_ = —0.5 result (the dashed curve). 




Figure 8: The Dalitz plots for the decay of ir^e^Ue with the square of matrix element 

summed over the electron polarizations (the left picture, po in arbitrary units) and with the 
spin polarization along the electron momentum (the right picture, in the same units as po). 
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Figure 9: The distribution of normal polarization on the Dahtz plot for the decays of K'^ — > 
n^ji^v^ (the left picture) and — >■ i^^e^Vf. (the right picture). The contours are shown with 
the arithmetic regressions for the polarization: the muon —0.0525 — j ■ 0.0525, the positron 
10-3- (-0.45 -j -0.60). 



4 Constraints on the leptoquark scales 

Under the determination of hadronic matrix elements of quark currents we derive the ratios 
of form factors due to the contribution of leptoquark interactions, 



fs V2 



f+{0) 16Gf\Vsu\ {rris - mu)mK A-Iq' 

h _ rriK 1 

/+(0) ~ 32G'f|K«| rriK* AIq' 



rriK-K 1 ^^^^ 



(55) 



where we have supposed the positive definiteness of Yukawa-constant products with respect 
to the mixing Vsu- Then we extract the values of leptoquark scales in the tensor part, 

Alq = 0.48+^17 TeV, (56) 

while the scalar form factor gives more stringent limit 

Alq = 3.4t~ TeV. (57) 

Thus, we deduce the 95%-confidence level 

Alq > 1.2 TeV. 

Let us compare the above restriction on the parameters of leptoquark interactions with the 
constraints following from other processes relevant to the effective vertices induced by diagrams 
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in Fig. Since the Yukawa constants are flavor dependent, the direct constraints can be 
obtained from the leptonic decays of kaon, viz., from both the electron and muon ones. In this 
way, the tensor interaction does not contribute, while the scalar one results in the multiplicative 
scaling of the decay amplitude. The factor has the form 

where we have neglected the masses of pion and w-quark, and mi denotes the mass of lepton. 
The leptonic modes are measured with the accuracy of branching ratios 

Be ~ 22' Bf, ~ 300' 

that can be used in order to restrict the scalar interactions induced by the leptoquarks. So, 
taking the ratio of branching ratios, which is independent of both the leptonic constant of 
kaon and the CKM element \Vsu\, we get the expression 

^ ^ fs rriK ^ ^ fs rriK 

_ ml /+(0) me _ ^ ^^^^ 5 /+(0) me 

/+(0) m, U{0) m, 

where we expand in small corrections following from the leptoquark interactions. Comparing 
with the experimental result 



Be 



= (2.44 ±0.11) ■ 10^ 

exp. 



we findQ 

Alq > 43 TeV. 

Thus, the measurements of semileptonic kaon decay provide us with the soft confirmation of 
constraints following from the leptonic decays, since the tensor and scalar effective vertices 
correlate in the leptoquark interactions. 



5 Discussion 

In this paper we have developed a model of meson dominance, which has allowed us to get 
quite an accurate description of hadronic form factors in the decay ir^l'^ui. In this 

way we have adjusted the model under the experimental data on the matrix element of vector 
quark-current and calculated the matrix element of tensor current induced by the leptoquark 

^As for some other restrictions see ref. [Q, where the bounds are very similar to those of obtained in the 
present paper. 
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interactions. The experimental data on the semileptonic decay of kaon allow us to extract the 
constraints on the contributions beyond the Standard Model, so that 

Alq > 1.2 TeV, 

where Alq represents the ratio of leptoquark mass to the square of Yukawa-like coupling. This 
limit softly confirms the bounds following from the leptonic decays of kaon. 

The authors thank the DELPHI and LHCb collaborations for a kind hospitality during 
their visit to CERN, where this paper was in progress. 

This work was in part supported by the Russian foundation for basic research, grants 
01-02-99315, 01-02-16585 and 00-15-96645, the Russian Ministry on the education, grant EOO- 
3.3-62. 



References 

[1] K. Horie et al. [KEK-E246 Collaboration], Phys. Lett. B 513, 311 (2001) 
[ |arXiv:hep-ex/0106006i ; 
A. S. Levchenko et al. [KEK-PS E246 Collaboration], |arXiv:hep-ex/0111048 



[2] I. V. Ajinenko et a/.. Preprint IHEP 2001-51, |arXiv:hep-ex/0lll>023| ; 
I. V. Ajinenko et al, Preprint IHEP 2002-6, |arXiv:hep-ph/ 020206 1| . 

[3] D. E. Groom et al. [Particle Data Group Collaboration], Eur. Phys. J. C 15, 1 (2000). 

[4] V. V. Braguta, A. A. Likhoded and A. E. Chalov, Phys. Rev. D 65, 054038 (2002) 
[ |arXiv:hep-ph/0106147| ; 

V. V. Braguta, A. A. Likhoded and A. E. Chalov, |arXiv: hep-ph / 1 05 1 11] ; 



R. N. Rogalev, Phys. Lett. B 521, 243 (2001) [|arXiv:hep-ph/01051871| ; 
Y. Kudenko, Part. Nucl. Lett. 106, 31 (2001). 

[5] K. Anikeev et al, |arXiv:hep-ph/0201071| ; 
V. V. Kiselev, |arXiv:hep-ph/0l09089| . 

[6] M. V. Chizhov, Phys. Lett. B 381, 359 (1996) ||arXiv:hep-ph/9511287 

[7] R. J. Tesarek [KTeV Collaboration] , |arXiv:hep-ex/9903069 



[8] W. BuchmuUer, R. Ruckl and D. Wyler, Phys. Lett. B 191, 442 (1987) [Erratum-ibid. B 
448, 320 (1987)]. 

[9] J. Kalinowski, R. Ruckl, H. Spiesberger and P. M. Zerwas, Z. Phys. C 74, 595 (1997) 
[|arXiv:hep-ph/9703288l . 



[10] G. Abbiendi et al [OPAL Collaboration], Phys. Lett. B 526, 233 (2002) 
[|arXiv:hep-ex/011202l . 



18 



[11] J. Gasser and H. Leutwyler, Nucl. Phys. B 250, 465 (1985), Nucl. Phys. B 250, 517 
(1985). 

[12] V. V. Kiselev, A. E. Kovalsky and A. I. Onishchenko, Phys. Rev. D 64, 054009 (2001) 
[|arXiv:hep-ph/0005020l . 

[13] S. Narison, |arXiv:hep-ph/ 0202200 . 

[14] S. S. Gershtein, A. A. Likhoded and A. I. Onishchenko, Phys. Kept. 320, 159 (1999). 



19 



